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ABSTRACT 


This  research  investigated  hlamentation  in,  and  modification  of,  dielectric  solids  by  multi- 
chromatic  ultrashort  laser  pulses.  Filamentation  of  ultrashort  laser  pulses  is  a  phenomenon 
with  a  growing  pool  of  applications  in  communications  and  national  defense.  Ultrafast  laser 
systems  can  also  affect  a  wide  variety  of  micrometer  and  nanometer  scale  modifications  to 
dielectrics  and  semiconductors.  This  research  was  a  theoretical  effort  to  develop  comprehen¬ 
sive  models  of  multi-chromatic  laser-material  interactions  in  bulk  solids  and  couple  them  to 
pulse  propagation  simulaisions.  The  research  investigated  techniques  to  control  filamentation 
and  the  onset  of  laser-induced  bulk  modification  with  multi-chromatic  pulses.  This  included 
controlling  modification  and  filamentation  in  solids  using  co-propagating  laser  pulses  of  dif¬ 
ferent  frequencies  as  well  as  strongly  chirped  pulses  and  spatially  chirped  annular  pulses. 
Emphasis  was  placed  on  nonlinear  optical  propagation  effects  since  they  determine  the  field 
distribution  at  the  laser  focus,  the  control  of  which  is  necessary  for  progress  in  laser  bulk 
micro-machining  towards  nano-machining.  Propagation  simulations  were  coupled  to  newly 
developed  extended  multi-rate  equation  methods  to  model  the  free-carrier  distribution  evolu¬ 
tion  in  the  energy  space  when  exposed  to  strongly  chirped  pulses.  Computationally  efficient 
3D  quantum  models  of  laser-solid  interactions  were  also  developed  for  coupling  with  the  ex¬ 
tended  multi-rate  equations  and  already  established  pulse  propagation  methods.  Coupling 
this  model  to  3D  pulse  propagation  and  comparing  the  results  to  a  collaborative  experiment 
conducted  at  AFRL  are  in  progress  at  the  time  of  this  report.  This  work  was  sponsored  by 
the  Air  Force  Office  of  Scientific  Research  (AFOSR)  within  the  Young  Investigator  Program. 


1  Summary 

1.1  Overview 

This  funding  fully  supported  a  postdoc  for  all  but  four  months  of  the  grant’s  duration,  as  well 
as  the  Pi’s  summer  salary  for  4.5  months.  In  addition  to  travel  and  publication  funds,  the 
grant  also  paid  for  one  travel  laptop,  one  visualization  desktop  unit,  and  two  computational 
servers.  The  research  supported  by  this  grant  is  disseminated  in  eight  publications,  one 
paper  currently  under  review,  two  papers  in  preparation,  and  six  conference  presentations. 
Each  of  these  works,  with  the  exception  of  the  collaborative  paper  under  review,  were  fully 
supported  by  this  funding  and  are  summarized  below. 

1.2  Publications 

•  J.  R.  Gulley,  J.  Liao,  and  T.  E.  Lanier,  “Plasma  generation  by  ultrashort  multi-chromatic 
pulses  during  nonlinear  propagation,”  Frontiers  in  Ultrafast  Optics:  Biomedical,  Scien¬ 
tific,  and  Industrial  Application  XIV ,  Proc.  SPIE  Int.  Soc.  Opt.  Eng.  8972,  89720T 
(2014). 
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J.  Liao  and  J.  R.  Gulley,  “Time- frequency  control  of  ultrafast  plasma  generation  in 
dielectrics,”  J.  Opt.  Soc.  Am.  B  31,  2973-2980  (2014). 

J.  R.  Gulley  and  T.  E.  Lanier, “Model  for  ultrashort  laser-pulse  induced  ionization 
dynamics  in  transparent  solids,”  Phys.  Rev.  B.  90,  155119  (2014). 

T.  E.  Lanier  and  J.  R.  Gulley,  “Annular  space-time  focusing  in  fused  silica,”  Frontiers 
in  Optics  2015,  OSA  Technical  Digest,  paper  FTu5E.6  (2015). 

J.  R.  Gulley  and  T.  E.  Lanier,  “Self-consistent  modeling  of  photoionization  and  the 
Kerr  effect  in  bulk  solids,”  Laser-Induced,  Damage  in  Optical  Materials:  2015,  Proc. 
SPIE  Int.  Soc.  Opt.  Eng.  7632,  96320X  (2015). 

T.  E.  Lanier  and  J.  R.  Gulley,  “Calculation  of  nonlinear  optical  damage  from  space- 
time-tailored  pulses  in  dielectrics,”  Laser-Induced  Damage  in  Optical  Materials:  2015, 
Proc.  SPIE  Int.  Soc.  Opt.  Eng.  9632,  96320Z  (2015). 

T.  E.  Lanier  and  J.  R.  Gulley,  “Nonlinear  space-time  focusing  and  filamentation  of 
annular  femtosecond  pulses  in  dielectrics,”  J.  Opt.  Soc.  Am.  B  33,  292-301  (2016). 

K.  Huthmacher,  A.  K.  Molberg,  B.  Rethfeld,  and  J.  R.  Gulley,  “A  C++ll/Fortran 
2003  interface  to  include  electron-electron  collisions  via  Monte  Carlo  in  Multiple  Rate 
Equation  simulations,”  J.  Comput.  Phys.  (under  review). 

J.  R.  Gulley  and  T.  E.  Lanier,  “Self-consistent  modeling  of  photoionization  and  non¬ 
linear  optics  in  cubic  dielectric  solids,”  (in  preparation). 

T.  E.  Lanier  and  J.  R.  Gulley,  “Two  pulse  multi-frequency  induced  modification  in 
optical  glasses,”  (in  preparation). 

Presentations 

J.  R.  Gulley,  J.  Liao,  and  T.  E.  Lanier,  “Plasma  generation  by  ultrashort  multi-chromatic 
pulses  during  nonlinear  propagation,”  Proc.  SPIE  8972,  Frontiers  in  Ultrafast  Optics: 
Biomedical,  Scientific,  and  Industrial  Application  XIV,  San  Francisco,  CA  (2014). 

J.  R.  Gulley,  “Self-consistent  modeling  of  photoionization  for  nonlinear  pulse  propa¬ 
gation  in  solids,”  Progress  in  Electromagnetics  Research  Symposium,  Prague,  Czech 
Republic  (2015). 

T.  E.  Lanier  and  J.  R.  Gulley,  “Calculation  of  nonlinear  optical  damage  from  space- 
time-tailored  pulses  in  dielectrics,”  Laser-Induced  Damage  in  Optical  Materials:  2015, 
Boulder,  CO  (2015). 

J.  R.  Gulley  and  T.  E.  Lanier,  “Self-consistent  modeling  of  photoionization  and  the 
Kerr  effect  in  bulk  solids,”  Laser-Induced  Damage  in  Optical  Materials:  2015,  Boulder, 
CO  (2015). 
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•  T.  E.  Lanier  and  J.  R.  Gulley,  “Annular  space-time  focusing  in  fused  silica,”  Frontiers 
in  Optics  2015,  San  Jose,  CA  (2015). 

•  J.  R.  Gulley,  “Self-consistent  modeling  of  photoionization  and  nonlinear  optics  in  di¬ 
electrics,”  High  Power  Laser  Ablation  /  Directed  Energy  Symposium,  Santa  Fe,  NM 
(2016). 

2  Introduction 

The  goal  of  this  project  was  to  model  nonlinear  ultrashort  laser  pulse  propagation  and  to  de¬ 
velop  techniques  that  make  novel  use  of  their  unique  multi-chromatic  properties.  Tradition¬ 
ally,  propagation  simulations  studying  ultrashort  pulse  hlamentation  and  bulk  modification 
include  an  internal  contradiction  for  the  case  of  solids:  they  combine  time-dependent  pulse 
amplitudes  with  the  monochromatic  approximation  underlying  the  electron-plasma  and  ion¬ 
ization  models.  Absence  of  an  adequate  model  that  accounts  for  the  multi-chromatic  nature 
of  ultrashort  laser  pulses  has  precluded  progress  in  the  applications  of  pulse  hlamentation,  as 
well  as  important  techniques  of  micromachining  and  nanofabrication  in  transparent  solids. 
This  project  set  out  to  couple  comprehensive  models  of  photoionization  and  multiphoton 
absorption  by  ultrashort  laser  pulses  with  a  newly  developed,  comprehensive  model  of  laser 
interactions  with  the  generated  conduction  electrons.  These  models  were  then  integrated 
into  the  modeling  of  pulse  propagation.  Additionally,  two  cases  of  interest  were  investigated 
for  controlling  hlamentation  and  bulk  modihcations  that  utilize  the  multi-chromatic  nature 
of  the  helds  (multi- frequency  multi-pulse  helds  and  chirped  pulses).  The  specihc  objectives 
of  the  project  were: 

1)  Develop  a  multi-chromatic  model  of  solid-state  plasma  dynamics  that  generates  a 
free-current  density  for  use  in  Maxwell’s  equations. 

2)  Integrate  this  model  into  a  numerical  code  to  simulate  ultrashort-pulse  propagation 
and  hlamentation  in  dielectric  solids. 

3)  Simulate  high-power  ultrafast  propagation  and  bulk  modihcation  by  multi-chromatic 
laser  pulses  in  solids  for  two  specihc  applications: 

(a)  strongly  chirped  pulses. 

(b)  two  co-propagating  pulses  of  different  frequencies. 

4)  Develop  non-time-averaged  multi-chromatic  quantum-mechanical  models  of  photoion¬ 
ization  and  the  Kerr  effect  in  dielectric  solids  by  ultrashort  laser  pulses. 

The  remainder  of  this  report  summarizes  the  accomplishments  and  progress  towards  these 
ends.  Where  the  details  of  research  hndings  are  available  in  published  form,  the  papers  are 
cited  in  the  section  titles  and  the  most  important  results  are  discussed.  Where  the  research 
is  not  yet  published,  more  detailed  descriptions  are  provided. 
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Figure  1:  On  the  left  side:  electron  energy  distributions  after  exposure  to  (a)  an  unchirped 
pulse,  (b)  a  negatively  chirped  pulse,  and  (c)  a  positively  chirped  pulse.  On  the  right  side: 
the  instantaneous  frequency  (d)  and  (e),  intensity  (f)  and  (g),  and  normalized  transmitted 
spectra  (h)  and  (i),  at  propagation  distances  of  z  =  0  and  20  m,  respectively.  The  negatively 
chirped  pulse  is  shown  as  the  blue  dashed  line,  the  unchirped  pulse  is  shown  as  the  black 
solid  line,  and  the  positively  pulse  is  shown  as  the  red  dotted  line. 
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3  Modeling  of  ionization  and  propagation  in  dielectrics 
by  strongly  chirped  pulses  [1-2] 

A  comprehensive  model  of  ultrafast  laser-induced  plasma  dynamics  intended  for  coupling 
with  pulse  propagation  simulations  in  transparent  solids  was  developed.  The  model  com¬ 
prised  two  sets  of  extended  multi-rate  equations  for  the  free-carrier  distribution  and  the  free- 
current  density.  The  changing  spectrum  of  a  propagating  ultrashort  laser  pulse  was  taken 
into  account  by  use  of  the  instantaneous  frequency  in  calculating  photoionization  and  single¬ 
photon  absorption  rates.  It  simultaneously  coupled  to  the  evolution  of  the  energy-resolved 
nonequilibrium  free-carrier  distribution,  which  was  examined  as  a  function  of  time  (see  left 
side  of  Fig.  1).  The  model  accounted  for  the  distribution  dependence  of  electron-phonon  and 
electron-electron  collision  times  in  the  distribution  evolution  as  well  as  free-carrier  absorp¬ 
tion,  each  using  relaxation  approximations. 

Our  simulations  solved  this  model  in  fused  silica  for  strongly  chirped  pulses  (shown 
Fig.  1).  The  results  suggest  that  strong  pulse  chirps,  such  as  those  occurring  naturally  dur¬ 
ing  nonlinear  propagation,  lead  to  ionization  dynamics  that  are  not  captured  by  the  stan¬ 
dard  monochromatic  treatment  of  laser-induced  plasma  formation.  The  order-of-magnitude 
chirp-dependent  differences  in  ionization  yield  and  the  electron  distribution  shapes  govern 
the  strength  and  character  of  the  laser-plasma  interactions.  The  proposed  model  therefore 
provides  a  framework  for  laser-material  interactions  that  is  practical  and  provides  better 
insight  into  the  physics  of  laser-induced  ionization  by  intense  ultrashort  laser  pulses.  Propa¬ 
gation  of  these  pulses  over  20  /i m  demonstrates  how  the  multi-chromatic  nature  of  oppositely 
chirped  pulses  alters  the  pulse  and  spectral  evolution  due  to  ionization.  Fully  3D  pulse  prop¬ 
agation  simulations  using  this  extended  multi-rate  equation  approach  also  demonstrates  a 
chirp  dependence  on  the  shape  and  localization  of  ionization  regions  in  the  material.  The  use 
of  the  instantaneous  frequency  (which  is  not  always  meaningful)  should  only  be  used  when 
a  central  frequency  can  still  be  well  defined.  The  intent  of  this  model  is  to  unify  it  with 
the  optical  Bloch  equation  methods  described  later  in  Sec.  7,  thus  providing  a  simultaneous 
quantum  mechanical  treatment  of  ionization  and  Kerr  effect  while  treating  the  many-body 
effects  using  the  relaxation  approximation. 


4  Ultrafast  ionization  in  dielectrics  with  stochastic  mod¬ 
eling  of  electron-electron  collisions  [3] 

In  addition  to  the  extended  multi-rate  equation  described  in  Sec.  3,  we  developed  a  split- 
step  numerical  method  for  calculating  ultrafast  free-electron  dynamics  in  a  collaboration 
with  the  group  of  Dr.  Baerbel  Rethfeld  at  the  University  of  Kaiserslautern  in  Germany. 
The  first  step  of  the  split-step  approach  solves  a  the  ordinary  deterministic  extended  multi¬ 
rate  equation  for  the  ionization,  electron-phonon  collisions,  and  single  photon  absorption 
by  free-carriers.  The  second  step  is  stochastic  and  models  electron-electron  collisions  by 
solving  the  Boltzmann  scattering  equation  using  Monte-Carlo  techniques,  thus  eliminating 
the  need  for  a  relaxation  treatment  of  electron-electron  collisions  which  are  the  dominant 
force  in  reshaping  the  non-equilibrium  distribution.  This  combination  of  deterministic  and 


5 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Energy  (eV) 

Figure  2:  Conduction  band  electron  distributions  in  fused  silica  as  functions  of  electron 
energy  at  various  times  for  400  nm  pulses  with  peak  intensities  of  4  x  1013  W/cm2  and 
durations  of  (a)  rw  =  15  fs,  (b)  rw  =  30  fs,  and  (c)  rw  =  60  fs. 
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stochastic  approaches  is  a  unique  and  efficient  method  of  calculating  the  nonlinear  dynamics 
of  3D  materials  exposed  to  high  intensity  ultrashort  pulses.  The  results  in  Fig.  2  from 
simulations  solving  the  split-step  model  in  fused  silica  demonstrate  that  electron-electron 
scattering  relaxes  the  non-equilibrium  electron  distribution  on  the  order  of  15  fs,  during 
the  exposure  to  ultrashort  pulses.  We  also  coupled  this  model  with  pulse  propagation,  and 
where  identical  pulses  and  material  parameters  are  assumed  the  results  agree  very  well  with 
those  of  the  extended  multi-rate  equation  using  the  relaxation  approximation,  which  is  much 
faster.  A  continued  collaboration  to  also  treat  the  important  impact  ionization  stochastically 
is  on-going. 


5  Simple  modeling  of  bulk  modification  by  multi- frequency, 
multi-pulse  fields  [4] 
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Figure  3:  Ionization  yields  generated  in  fused  silica  by  50  fs  -  250  fs  double-pulse  trains  as 
functions  of  time  (given  in  the  reference  frame  of  the  250  fs  pulse)  and  delay  between  the 
two  pulses.  The  50  fs  and  250  fs  pulses  have  wavelengths  as  labeled  and  respective  peak 
intensities  of  1  x  1012  W/cm2  and  4  x  1012  W/cm2. 


In  the  recent  ultrashort  laser-induced  micromachining  literature  there  has  been  increased 
interest  in  experiments  using  multi-color,  two-pulse  fields.  Typically,  this  involves  using  a 
high  frequency,  low  energy  pulse  to  drive  photoionization  followed  by  a  low  frequency,  high 
energy  pulse  to  drive  avalanching.  To  date,  the  experimental  papers  using  this  technique 
attempt  to  verify  this  approach  by  modeling  the  ionization  with  oversimplified  ID  prop¬ 
agation  models  combined  with  one  or  two  rate  equations  to  calculate  the  ionization.  To 
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facilitate  a  fast  and  effective  modeling  for  time-frequency  control  of  plasma  generation  for 
fully  3D  propagation  simulations,  we  developed  a  modified  ordinary  multi-rate  equation  to 
calculate  ionization  by  double-pulses  of  different  frequencies.  This  model  also  discretizes  the 
conduction  electron  distribution  in  energy  space,  but  with  levels  equal  to  the  smallest  photon 
energy  of  the  low-frequency  pulse.  The  approach  is  comparable  in  computational  speed  to 
the  ubiquitous  single  rate  equation  for  ionization  in  the  filamentation  literature.  We  used 
this  approach  to  calculate  the  optimal  experimental  delay  between  pulses  of  particular  fre¬ 
quencies,  energies  and  pulsewidths  (see  Fig.  3).  In  practice,  the  first  of  these  pulses  is  usually 
a  harmonic  of  the  second  pulse  at  the  fundamental  frequency.  However,  the  Fig.  3  results 
on  the  right  side  show  that  non-harmonic  double-pulse  trains  can  control  plasma  generation 
with  a  minimal  fluence  while  still  using  common  laser  frequencies. 


6  Control  of  filamentation  and  bulk  modification  by 
spatio-temporally  chirped  pulses  [5-7] 
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Figure  4:  On  the  left:  a  schematic  of  the  simulated  SSTF  system  comprising  two  concentric 
gratings  Gi  and  G2  with  ring  spacing  d  =  10  //in  and  a  lens  of  focal  length  /  =  200  mm.  On 
the  right:  ionization  tracks  for  0.5-//J  (a)  Gaussian  beam-pulse  with  initial  an  pulsewidth  of 
40  fs  and  (b)  SSTF  annular  pulse  of  approximately  40  fs  pulsewidth  at  the  linear  focus. 


Filamentation  during  simultaneous  space-time  focusing  in  bulk  fused  silica  was  investi¬ 
gated  numerically  for  spatially  chirped  annular  pulses.  We  modeled  the  use  of  a  pair  of 
concentric  gratings  to  transform  a  common  femtosecond  laser  pulse  with  Gaussian  spatial 
profile  into  a  radially  chirped,  annular  beam  shape  that  is  focused  into  the  bulk  of  silica 
by  a  lens.  The  system  schematic  for  this  setup  is  shown  on  the  left  side  of  Fig.  [4],  Pulse 
propagation  simulations  compared  the  ionization  yields  of  these  annular  space-time  focus¬ 
ing  waveforms  as  functions  of  beam  radius  and  propagation  distance,  and  are  compared  to 
those  of  ordinary  gaussian  pulses  of  the  same  numerical  aperture,  shown  on  the  right  side 
of  Fig.  [4],  The  results  show  that  slightly  higher  ionization  yields,  and  greater  radial  lo¬ 
calization  of  the  peak  structure  of  these  yields  on  the  micro-nanometer  scale,  are  possible 
with  the  annular  space-time  focusing  configuration.  This  is  significant  since  the  regions  of 
peak  ionization  yield  often  correlate  with  regions  of  laser-induced  material  modification  in 
multi-shot  experiments. 
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Figure  5:  Left  side  shows  SSTF  in  vacuum  of  annular  pulses  arising  from  radial  chirp, 
(a)  Tilted  rco-spectrum  just  after  G2.  (b)-(c)  Normalized  intensities  at  various  ^-positions 
(superposed).  In  each  case,  propagation  distances  as  measured  from  the  /  =  200  mm  lens  are 
0  mm  (top),  100  mm  (middle),  and  195  mm  (bottom),  (b)  The  initially  bandwidth- limited 
40  fs  pulse,  (c)  The  positively  pre-chirped  pulse.  Dashed  lines  are  a  guide  for  the  eye.  Right 
side  shows  the  annular  SSFT  pulse  just  beyond  the  nonlinear  focus  for  a  0.7-//J  pulse  with 
optimal  SSTF  pre-chirp,  (d)  The  log  of  the  intensity  (W/m2)  in  rf-space.  (e)  Normalized 
intensity  (inset),  and  the  log  of  the  normalized  intensity  (6  decades)  in  fc^cu-space. 


It  is  typically  assumed  that  in  simultaneous  space-time  focusing  (SSTF)  configurations 
the  nonlinear  propagation  effects  are  approximately  negligible.  Our  results  confirm  this  out¬ 
side  of  the  nonlinear  focal  region,  but  inside  the  nonlinear  focus  the  filamentation  dynamics 
during  the  SSTF  propagation  play  a  non-negligible  role  in  the  pulse  and  material  evolution. 
To  acquire  the  optimal  SSTF  condition  of  a  simultaneous  minimum  in  both  beam  and  pulse 
widths,  one  must  vary  the  energy  and/or  time-chirp  of  the  incident  pulse  (see  left  side  of 
Fig.  [5]).  As  the  annular  pulse  propagates  in  the  bulk  silica  we  capture  the  field  dynamics 
and  material  response  that  yields  material  modification  and  damage  in  the  nonlinear  focus. 
The  propagation  results  show  rich  nonlinear  dynamics  typical  of  filamenting  pulses,  see  right 
side  of  Fig.  [5].  Thus,  while  higher  ionization  yields  and  greater  spatial  localization  along 
the  beam  radius  of  these  yields  were  obtained  with  annular  SSTF  pulses,  such  pulses  were 
not  effective  at  restricting  the  filamentation  dynamics  to  the  same  length  scales  along  the 
propagation  axis. 
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7  Quantum  modeling  of  photoionization  and  nonlinear 
optics  in  cubic  solids  [8] 

In  this  section  we  describe  results  for  modeling  the  laser-induced  electron  dynamics  for  a 
simple  cubic  3D  solid  using  a  reduced  set  of  Optical  Bloch  equations.  It  was  important  to 
move  beyond  the  two-band  approach  originally  proposed  in  the  grant  in  order  to  calculate  the 
Kerr  effect  and  photoionization  self-consistently.  This  approach  calculates  the  contributions 
of  photoionization,  the  Kerr  effect,  and  current  from  free  carriers  in  a  self-consistent  way 
that  makes  no  assumptions  about  the  pulse  spectrum.  In  the  perturbative  nonlinear  optical 
regime,  the  resulting  polarization  from  this  approach  accurately  reproduces  the  linear  and 
nonlinear  optical  properties  of  wide-band  gap  cubic  solids  with  corresponding  multi-photon 
ionization.  At  higher  field  strengths  corresponding  to  the  non-perturbative  nonlinear  optical 
(tunneling)  regime  the  polarization  includes  high-order  harmonic  generation  due  both  to 
high  order  Kerr  effects  and  Bragg  reflection  in  the  crystal.  This  approach  calculates  all 
orders  of  nonlinear  optical  polarizations  resulting  from  electronic  transitions  from  valance  to 
conduction  bands.  The  case  of  photoionization  by  multiple  laser  frequencies  was  examined 
in  particular.  Also,  photoionization  and  the  electronic  Kerr  effect  are  shown  to  be  strongly 
coupled  for  ultrashort  pulses  on  the  order  of  20  fs  or  less,  contrary  to  the  predictions  of  more 
traditional  approaches.  This  is  on  the  same  time  scale  as  the  electron-electron  collisions 
discussed  in  Sec.  4,  and  the  intent  in  developing  this  approach  is  to  merge  it  with  the 
extended  multi-rate  equation  for  a  comprehensive  model  of  ionization  and  nonlinear  optics. 
These  results  are  still  in  preparation  for  publishing  at  the  time  of  this  report,  so  we  provide 
theoretical  detail  sufficient  to  reproduce  the  presented  results. 


7.1  Multi-Band  Laser-Solid  Dynamics 

The  optical  transitions  in  a  multi-band  3D  quasi-momentum  space  k  can  be  calculated  by  the 
Semiconductor  Bloch  Equations  (SBEs).  On  the  femtosecond  time  scale  the  optical  effects 
of  many-body  Coulomb  processes  can  be  modeled  by  a  phenomenological  dephasing  rate  7, 
effectively  reducing  the  Semiconductor  Bloch  Equations  to  a  set  of  Optical-Bloch  Equations. 
The  time  evolution  of  the  occupation  numbers  for  electrons  (/e„,k)  and  holes  (fhn, k),  and 
the  corresponding  microscopic  polarizations  were  generalized  for  multiband  systems  in  the 
supplementary  material  of  [O.  Schuber,  et  al,  Nat.  Photon.  8,  119  (2014)]: 
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Here,  e^m5k  and  e6n?k  are  the  energies  of  the  holes  and  electrons  in  the  mth  and  nth  bands, 
respectively.  The  terms  with  dipole  matrix  elements  dXiVj^  cause  the  interband  transitions 
while  the  terms  with  zV k  excite  intraband  activity.  The  laser  field  E(t )  is  assumed  to  be 
linearly  polarized.  This  model  allows  the  inclusion  of  an  arbitrary  number  of  electron  and 
hole  bands,  but  the  band  structure  must  still  be  specified. 

7.2  Band  Structure 

Our  model  is  developed  for  laser-material  interactions  in  a  simple  cubic  solid.  For  the  sake  of 
simplicity  while  still  retaining  desirable  band  structure  properties  (i.e.  non-parabolic  shape, 
periodicity  in  quasi-momentum  space,  etc.)  we  use  analytic  formalisms  to  approximate  the 
full  3D  energy  dispersion  relations  to  ensure  these  properties.  One  possible  model  is  the 
energy  relation  for  a  simple  cubic  medium  in  the  nearest  neighbor  tight  binding  limit: 

£n,  k  =  £n,r  +  An  (1  -  COS  (kid))  ,  (1) 
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where  en  r  is  the  electron  energy  of  the  nth  band  at  the  Gamma  point.  The  coefficients 
An  =  h2/(mnd2),  d  is  the  lattice  constant,  and  mn  is  the  diagonal  component  of  the  electron 
effective  mass  tensor  at  k  =  0  for  the  nth  band. 

Equation  1  is  an  instructive  dispersion  relation  that  we  used  to  model  a  hypothetical 
9  eV  band  gap  crystal.  However,  to  simulate  a  more  realistic  solid  we  chose  the  simple 
cubic  crystal  Csl.  The  nearest  neighbor  tight  binding  relation  in  Eq.  1  can  be  fitted  to  well 
approximate  the  valence  bands  in  Csl  when  compared  to  theoretical  calculations.  Eq.  1 
also  fits  well  with  the  Csl  conductions  bands  between  the  T-point  and  the  X-point,  but  not 
between  the  T-point  and  the  M-points  or  R-points.  A  better  fit  (although  still  not  perfect) 
for  the  conduction  bands  is  provided  by  the  relation: 

£n, k  =  £n,r  +  An  -  An  cos (kxd)  cos (kyd)  cos (kzd),  (2) 

which  has  been  used  to  approximate  the  transition  energies  in  Alkalide  Halides  in  the  liter¬ 
ature  before  [V.  E.  Gruzdev,  Phys.  Rev.  B  75,  205106  (2007)].  We  therefore  use  Eq.  2  to 
approximate  the  structure  of  the  Csl  bands.  All  energies  at  the  T-points  are  taken  from  the 
band  calculations  of  [Y.  Onodera,  J.  Phys.  Soc.  Jap.  25,  469  (1968)],  and  the  values  for  An 
are  approximated  from  the  calculated  band  structures  by  An  =  (en,x  ~  Gi,r)/ 2-  From  this 
value  we  extract  the  approximate  effective  mass  at  each  T-point,  mn. 

7.3  Selection  of  dipole- matrix  elements 

Solving  the  SBEs  with  the  desired  band  structure  requires  dipole-matrix  elements  (%k, 
which  can  be  appreciably  affected  by  the  chosen  band  structure.  In  our  calculations  we  vary 
the  dipole  coupling  constant  with  momentum  according  to  the  well  know  relation  dcv^  = 
dcv, o(c  r~£v  r)/(£c,k— £v,k),  but  this  still  requires  selecting  a  value  for  dcv  (}.  While  accounting 
first  for  necessary  selection  rules,  the  remaining  values  of  dCV: o  are  estimated  by  k  •  p  theory 
and  then  scaled  to  return  the  linear  susceptibility  as  given  by  the  sum  over  k-states 

(i) \  ■>  f  4  |dct,  k|  £cv, k  d?k 

U  “Wbz  (27r)3e0  (e«,.k  -  ^  -  *' 7) (c*,,k  T  huj  +  hy) 

The  remaining  dipole  matrix  elements  are  likewise  chosen  to  return  a  reasonable  value  for 
the  third-order  susceptibility  for  long  wavelengths  10-18cm2/V2. 

7.4  Reduction  to  effective  2D  Calculation 

The  SBEs  for  a  3D  solid  can  be  solved  directly  in  the  3D  quasi-momentum  space  of  the 
crystal,  but  this  is  computationally  impractical  for  long  distance  propagation  simulations.  If 
the  crystal  structure  is  centro-symmetric  and  the  laser  field  is  linearly  polarized  along  one 
of  the  crystal  axes  (assumed  to  be  the  ^-direction),  this  procedure  leads  to  many  identical 
calculations  that  provide  no  new  information.  This  is  because  the  intrinsic  inversion  sym¬ 
metry  ensures  that  ,/n ( kx ,  ky ,  kz )  f n (  kX)  ky ,  kz )  and  P\ j ( />y ,  ky .  kz )  P\j (  />y ,  ,  kz ) . 

One  may  not  assume  this  inversion  symmetry  in  all  three  directions  since  the  presence  of  the 
laser  field  breaks  the  inversion  symmetry  in  the  direction  of  the  field. 
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If  the  uniqueness  of  each  calculation  in  a  kx-ky  plane  at  a  particular  kz  is  determined  by 
its  associated  electron  (hole)  energy  e(k),  then  one  need  only  solve  each  unique  SBE  and  have 
knowledge  of  the  effective  density  of  states  for  that  solution.  To  accomplish  this  reduction  we 
transformed  the  (kx,  ky,  kz)  coordinate  system  to  an  effective  2D  ( Qxy ,  kz )  system.  Here  the 
dimensionless  coordinate  Qxy  is  a  dimensionless  coordinate  that  ranges  from  —  1  <  Qxy  <  1. 
The  energies  for  the  two  band  structures  are  then  given  by 

£n,k  £n,r  T  3An  A n  cos(kzd)  T  2/kt^^y,  (3) 


£n, k  &n, r  T  A n  A n  COS (kzd)Qxy , 


(4) 


In  either  case,  the  effective  density  of  states  under  this  new  coordinate  system  is  given  by 
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1  -Q2 


xy 


which  includes  electron  spin  degeneracy.  The  function  K  [X ]  is  the  complete  elliptical  integral 
of  the  first  kind.  With  this  simplification  the  SBEs  for  the  full  3D  momentum  space  can 
be  solved  in  real  time  on  the  order  of  minutes  (attempts  to  make  futher  simplification  to 
an  effective  lD-array  calculation  are  on-going  and  would  reduce  the  calculation  time  to  the 
order  of  1  s).  Under  this  new  coordinate  system,  the  total  ionization  yield  Ne ,  free  current 
density  -If ,  and  total  material  polarization  P  are  calculated  by: 


Ne  =  Y  I  fen(Q  x,yi  kz  )o{Qx,yi  kz)dQxydkz, 
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where  jn(k)  =  (|e|  /^)Vken(fc). 


7.5  Results  for  a  three-band  9-eV  simple  cubic  dielectric 

The  SBEs  were  solved  for  a  hypothetical  9-eV  simple  cubic  solid  in  the  nearest  neighbor  tight 
binding  limit.  The  three  bands  comprised  a  flat  p-type  valance  band,  one  s-type  conduction 
band,  and  one  higher  f-type  conduction  band.  Parameters  were  chosen  such  that  the  linear 
and  nonlinear  refractive  indices  were  n0  ~  1.5  and  n2  ~  3  x  10“ 16  W/cm2,  respectively.  For 
exposure  to  an  800  nm,  20  fs  pulse  with  a  peak  intensity  of  10  TW/cm2,  Fig.  6  shows  the 
final  occupation  numbers  in  the  first  conduction  band  as  a  function  of  Qxy  and  kz.  Here  one 
can  see  the  evidence  of  multi-photon  absorption  in  the  various  rungs  of  the  final  electron 
distribution.  We  note  that  there  were  also  electrons  in  the  second  conduction  band,  albeit 
with  occupation  numbers  an  order  of  magnitude  or  less  than  those  of  the  first  conduction 
band.  The  corresponding  nonlinear  polarization  is  shown  in  Fig.  7.  Also  shown  in  Fig.  7  is 
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Figure  6:  The  occupation  probability  fei  as  a  function  of  the  transverse  coordinate  Qxy  and 
the  momentum  kz  after  exposure  to  a  20  fs,  800  nm  laser  pulse  polarized  in  the  ^-direction. 


Figure  7:  The  calculated  nonlinear  polarization  (solid-black  line)  and  the  traditional 
monochromatic  Kerr-type  polarization  (dashed-red  line)  during  exposure  to  a  20  fs,  800 
nm  laser  pulse. 
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Figure  8:  The  calculated  free  current  density  (a)  and  free  carrier  polarization  (b)  as  calculated 
by  the  SBEs  (solid-black  lines)  and  the  classical  free  carrier  current  (dashed-red  line)  during 
exposure  to  a  20  fs,  800  nm  laser  pulse. 


the  calculation  of  the  traditional  Kerr  effect  Tkerr(^)  =  eo x^E3(t).  The  agreement  between 
the  two  models  is  initially  good.  As  ionization  increases,  however,  higher  harmonics  emerge 
in  the  full  calculation. 

Figure  8(a)  shows  the  full  current  density  calculation  along  with  the  classical  free  current 
density  calculated  by  Jc(t)  =  (e2 /mr)  p  E(t)Nk(t)dt,  where  Nk(t)  is  the  ionization  yield 
calculated  by  the  Keldysh  model.  Figure  8(b)  show  the  time  integrated  values  of  both 
quantities  in  Figure  8(a)  so  that  they  appear  as  polarization  contributions.  Initially  the  full 
free  polarization  makes  a  contribution  that  is  indistinguishable  from  the  ordinary  Kerr  effect. 
Only  after  ionization  occurs  does  the  full  free  current  match  the  form  of  the  classical  plasma 
current,  which  in  turn  leads  to  plasma  defocusing  of  a  laser  beam,  thus  demonstrating  that 
the  separation  of  Kerr  and  ionization  effects  is  artificial  for  intensities  typical  of  filamentation 
and  laser- induced  damage  in  solids. 

7.6  Results  for  seven-band  Csl  exposed  to  multi- harmonic  fields 

For  a  7-band  Csl  simple  cubic  crystal  we  solved  the  SBEs  for  20  fs  pulses  using  Eq.  2  for  the 
band  structure.  The  three  top  valence  bands,  the  three  lowest  conduction  bands,  and  the 
lowest  f-type  conduction  band  were  chosen  to  represent  the  band  structure.  For  the  800  nm 
pulse  in  isolation,  the  post-exposure  occupation  numbers  of  the  three  lowest  conduction 
bands  are  shown  in  Figure  9(a-c).  The  multi-photon  ionization  structure  of  these  plots  are 
difficult  to  interpret  until  one  compares  it  to  the  conduction  band  energies  as  a  function  of 
the  same  coordinates,  see  Figure  9(d). 
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Figure  9:  The  occupation  probabilities  (a)  fei,  (b)  /e2,  and  (c)  fea  as  a  function  of  the 
transverse  energy  Qxy  and  the  field-polarized  momentum  kz  after  exposure  to  a  20  fs,  800 
nm  laser  pulse,  (d)  The  conduction  band  energy  as  a  function  of  Qxy  and  kz . 


We  also  solved  the  SBEs  for  other  20  fs  single  pulses  and  multi-color  pulse  combinations. 
One  such  case  was  for  an  800  nm  pulse  of  peak  intensity  5  TW/cm2,  a  267  nm  pulse  of 
peak  intensity  0.05  TW/cm2,  and  both  pulses  present  simultaneously.  All  the  pulses  are 
assumed  to  be  Gaussian  and  centered  about  t  =  0.  Figure  10  shows  the  conduction  electron 
populations  as  functions  of  time  for  each  of  these  cases,  the  final  population  constituting 
the  net  yield.  The  267  nm  pulse  in  isolation  had  both  the  lowest  intensity  and  the  lowest 
yield,  with  the  yield  of  the  800  nm  pulse  being  an  order  of  magnitude  above  it.  However, 
the  photoionization  yield  for  the  combined  field  increases  by  an  order  of  magnitude,  which 
is  purely  a  multi-color  photoionization  effect  which  could  be  detectable  in  experiment.  Non¬ 
linear  propagation  simulations  with  such  field  combinations  is  challenging  because  the  band 
gap  of  Csl  (6.2  eV)  straddles  a  two-photon  absorption  for  the  267  nm  pulse,  and  using  a 
constant  nonlinear  refractive  index  is  not  valid  for  that  field  combination.  Figure  10(c)  shows 
the  nonlinear  polarization  for  the  combined  800  nm  and  267  nm  field  in  the  spectral  domain. 
Here  we  note  the  presence  of  both  even  and  odd  order  harmonics  of  the  fundamental,  which 
would  not  be  predicted  by  a  traditional  Kerr  effect  calculation  for  this  combind  field.  This 
is  due  primarily  to  the  presence  of  multiple  bands  and  would  also  not  be  captured  by  a 
two-band  model. 
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Figure  10:  (a)  The  independent  800  nm  (red)  and  267  nm  (purple)  pulse  intensities,  (b)  The 
total  ionization  yields  for  each  individual  pulse  and  for  the  combined  fields  (black),  (c)  The 
nonlinear  polarization  spectrum  for  the  combined  field  in  Csl. 


8  Conclusions 

This  research  investigated  filamentation  in,  and  modification  of,  dielectric  solids  by  multi- 
chromatic  ultrashort  laser  pulses  through  modeling  and  simulation.  We  simulated  techniques 
to  control  filamentation  and  the  onset  of  laser-induced  bulk  damage  with  multi-chromatic 
pulses  using  strongly  chirped  pulses  and  spatially  chirped  annular  pulses.  We  developed  ex¬ 
tended  multi-rate  equation  models  for  free-carrier  dynamics  of  sufficient  simplicity  to  couple 
with  fully  3D  pulse  propagation  simulations  while  accounting  for  many-body  collisions  using 
relaxation  approximations  and  stochastic  solving  of  a  Boltamann  scattering  equation.  We 
also  developed  a  computationally  efficient  3D  quantum  calculation  of  laser-solid  coupling  by 
solving  a  reduced  set  of  Semiconductor  Bloch  Equations  to  obtain  the  source  terms  for  prop¬ 
agation  models  and  examined  the  solutions  of  this  model  for  multi-chromatic  pulse  systems. 
This  work  is  disseminated  in  eight  papers  (one  under  review),  six  conference  presentations, 
and  two  papers  in  preparation  at  the  time  of  this  report. 


9  Bibliography 

1.  J.  R.  Gulley,  J.  Liao,  and  T.  E.  Lanier,  “Plasma  generation  by  ultrashort  multi-chromatic 
pulses  during  nonlinear  propagation,”  Frontiers  in  Ultrafast  Optics:  Biomedical,  Scien¬ 
tific,  and  Industrial  Application  XIV ,  Proc.  SPIE  Int.  Soc.  Opt.  Eng.  8972,  89720T 
(2014). 


17 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


2.  J.  R.  Gulley  and  T.  E.  Lanier, “Model  for  ultrashort  laser-pulse  induced  ionization 
dynamics  in  transparent  solids,”  Phys.  Rev.  B.  90,  155119  (2014). 

3.  K.  Huthmacher,  A.  K.  Molberg,  B.  Rethfeld,  and  J.  R.  Gulley,  “A  C++ 11 /Fortran 
2003  interface  to  include  electron-electron  collisions  via  Monte  Carlo  in  Multiple  Rate 
Equation  simulations,”  J.  Comput.  Phys.  (under  review). 

4.  J.  Liao  and  J.  R.  Gulley,  “Time- frequency  control  of  ultrafast  plasma  generation  in 
dielectrics,”  J.  Opt.  Soc.  Am.  B  31,  2973-2980  (2014). 

5.  T.  E.  Lanier  and  J.  R.  Gulley,  “Annular  space-time  focusing  in  fused  silica,”  Frontiers 
in  Optics  2015,  OSA  Technical  Digest,  paper  FTu5E.6  (2015). 

6.  T.  E.  Lanier  and  J.  R.  Gulley,  “Calculation  of  nonlinear  optical  damage  from  space- 
time-tailored  pulses  in  dielectrics,”  Laser-Induced  Damage  in  Optical  Materials:  2015, 
Proc.  SPIE  Int.  Soc.  Opt.  Eng.  9632,  96320Z  (2015). 

7.  T.  E.  Lanier  and  J.  R.  Gulley,  “Nonlinear  space-time  focusing  and  filamentation  of 
annular  femtosecond  pulses  in  dielectrics,”  J.  Opt.  Soc.  Am.  B  33,  292-301  (2016). 

8.  J.  R.  Gulley  and  T.  E.  Lanier,  “Self-consistent  modeling  of  photoionization  and  the 
Kerr  effect  in  bulk  solids,”  Laser-Induced  Damage  in  Optical  Materials:  2015,  Proc. 
SPIE  Int.  Soc.  Opt.  Eng.  7632,  96320X  (2015). 


18 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Response  ID:6185  Data 


1. 

1 .  Report  Type 
Final  Report 
Primary  Contact  E-mail 

Contact  email  if  there  is  a  problem  with  the  report. 

jgulley@kennesaw.edu 

Primary  Contact  Phone  Number 

Contact  phone  number  if  there  is  a  problem  with  the  report 

470-578-2933 

Organization  /  Institution  name 

Kennesaw  State  University 

Grant/Contract  Title 

The  full  title  of  the  funded  effort. 

(YIP)  Multi-Chromatic  Ultrashort  Pulse  Filamentation  and  Bulk  Modification  in  Dielectrics 

Grant/Contract  Number 

AFOSR  assigned  control  number.  It  must  begin  with  MFA9550M  or  MF49620M  or  "FA2386". 

FA9550-1 3-1 -0069 

Principal  Investigator  Name 

The  full  name  of  the  principal  investigator  on  the  grant  or  contract. 

Jeremy  R.  Gulley 
Program  Manager 

The  AFOSR  Program  Manager  currently  assigned  to  the  award 
Dr.  Arje  Nachman 
Reporting  Period  Start  Date 

02/01/2013 

Reporting  Period  End  Date 

01/31/2016 

Abstract 

This  research  investigated  filamentation  in,  and  modification  of,  dielectric  solids  by  multi-chromatic 
ultrashort  laser  pulses.  It  was  a  theoretical  effort  to  develop  models  of  multi-chromatic  laser-material 
interactions  in  bulk  solids  and  couple  these  to  pulse  propagation  simulations.  The  research  developed 
techniques  to  control  filamentation  and  the  onset  of  laser-induced  bulk  damage  with  multi-chromatic  fields. 
This  included  controlling  modification  and  filamentation  in  solids  using  co-propagating  laser  pulses  of 
different  frequencies  as  well  as  strongly  chirped  pulses  and  spatially  chirped  annular  pulses.  Emphasis 
was  placed  on  nonlinear  optical  propagation  effects  since  they  determine  the  field  distribution  at  the  laser 
focus,  the  control  of  which  is  necessary  for  progress  in  laser  bulk  micro-machining  towards  nano¬ 
machining.  Computationally  efficient  3D  quantum  models  of  laser-induced  photoionization  and  the 
electronic  Kerr  effect  were  developed  for  coupling  with  the  already  established  pulse  propagation  methods. 

Distribution  Statement 

This  is  block  12  on  the  SF298  form. 

Distribution  A  -  Approved  for  Public  Release 

Explanation  for  Distribution  Statement 

If  this  is  not  approved  for  public  release,  please  provide  a  short  explanation.  E.g.,  contains  proprietary  information. 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


SF298  Form 

Please  attach  your  SF298  form.  A  blank  SF298  can  be  found  here.  Please  do  not  password  protect  or  secure  the  PDF 
The  maximum  file  size  for  an  SF298  is  50MB. 

S98.pdf 

Upload  the  Report  Document.  File  must  be  a  PDF.  Please  do  not  password  protect  or  secure  the  PDF .  The 
maximum  file  size  for  the  Report  Document  is  50MB. 

final-report.pdf 

Upload  a  Report  Document,  if  any.  The  maximum  file  size  for  the  Report  Document  is  50MB. 

Archival  Publications  (published)  during  reporting  period: 

J.R.  Gulley,  J.  Liao,  and  T.E.  Lanier,"  Plasma  generation  by  ultrashort  multi-chromatic  pulses  during 
nonlinear  propagation,"  Frontiers  in  Ultrafast  Optics:  Biomedical,  Scien-  tific,  and  Industrial  Application  XIV, 
Proc.  SPIE  Int.  Soc.  Opt.  Eng.  8972,  89720T  (2014). 

J.  R.  Gulley  and  T.  E.  Lanier, "Model  for  ultrashort  laser-pulse  induced  ionization  dynamics  in  transparent 
solids,"  Phys.  Rev.  B.  90, 1 551 1 9  (2014). 

J.  Liao  and  J.  R.  Gulley,  "Time-frequency  control  of  ultrafast  plasma  generation  in  dielectrics,"  J.  Opt.  Soc. 

Am.  B  31, 2973-2980  (2014). 

T.  E.  Lanier  and  J.  R.  Gulley,  "Annular  space-time  focusing  in  fused  silica,"  Frontiers  in  Optics  2015,  OSA 
Technical  Digest,  paper  FTu5E.6  (2015). 

T.  E.  Lanier  and  J.  R.  Gulley,  "Calculation  of  nonlinear  optical  damage  from  space-  time-tailored  pulses  in 
dielectrics,"  Laser-Induced  Damage  in  Optical  Materials:  2015,  Proc.  SPIE  Int.  Soc.  Opt.  Eng.  9632,  96320Z 
(2015). 

J.  R.  Gulley  and  T.  E.  Lanier,  "Self-consistent  modeling  of  photoionization  and  the  Kerr  effect  in  bulk  solids," 
Laser-Induced  Damage  in  Optical  Materials:  2015,  Proc.  SPIE  Int.  Soc.  Opt.  Eng.  7632,  96320X  (2015). 

T.  E.  Lanier  and  J.  R.  Gulley,  "Nonlinear  space-time  focusing  and  filamentation  of  annular  femtosecond 
pulses  in  dielectrics,"  J.  Opt.  Soc.  Am.  B  33,  292-301  (2016). 

K.  Pluthmacher,  A.  K.  Molberg,  B.  Rethfeld,  and  J.  R.  Gulley,  "A  C++1 1 /Fortran  2003  interface  to  include 
electron-electron  collisions  via  Monte  Carlo  in  Multiple  Rate  Equation  simulations,"  J.  Comput.  Phys. 

(under  review). 

Changes  in  research  objectives  (if  any): 

Change  in  AFOSR  Program  Manager,  if  any: 

Extensions  granted  or  milestones  slipped,  if  any: 

AFOSR  LRIR  Number 
LRIR  Title 
Reporting  Period 
Laboratory  Task  Manager 
Program  Officer 
Research  Objectives 
Technical  Summary 

Funding  Summary  by  Cost  Category  (by  FY,  $K) 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Starting  FY 

FY+1 

FY+2 

Salary 

Equipment/Facilities 

Supplies 

Total 

_ 

_ 

Report  Document 

Report  Document  -  Text  Analysis 

Report  Document  -  Text  Analysis 

Appendix  Documents 

2.  Thank  You 

E-mail  user 

Apr  30, 2016  09:32:00  Success:  Email  Sent  to:  jgulley@kennesaw.edu 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


